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Abstract. - The ground state properties and low- lying excitations of a (quasi) one-dimensional 
system of longitudinally confined interacting bosons are studied. This is achieved by extending 
Haldane's harmonic-fluid description to open boundary conditions. The boson density, one- 
particle density matrix, and momentum distribution are obtained accounting for finite-size and 
boundary effects. Friedel oscillations are found in the density. Finite-size scaling of the momen- 
tum distribution at zero momentum is proposed as a method to obtain from the experiment 
the exponent that governs phase correlations. The strong correlations between bosons induced 
by reduced dimensionality and interactions are displayed by a Bijl-Jastrow wave function for 
the ground state, which is also derived. 



The names of Tomonaga and Luttinger are commonly associated with interacting elec- 
tron systems in one-dimension (Id). The so-called Tomonaga-Luttinger liquids (TLL's) have 
the remarkable property that their low-energy spectrum is completely exhausted by gapless 
collective excitations, i.e. there are no quasiparticles resembling the constituent electrons. 
However, as emphasized by Haldane the class of TLL's is broader than the electron sys- 

tems that have attracted much attention, especially in recent years. In this paper we shall be 
concerned with Id systems whose constituent particles are not electrons, but bosonic atoms. 
A few experimental systems which fit well these requirements are already available: Atomic 
vapors confined in highly anisotropic traps 1^ and the axial phase of ^He absorbed in narrow 
pores Q or nanotubes are good examples. The possibility of TLL behavior in trapped 
vapors has already been discussed in the literature |^,D. Monien et al. Q have pointed out 
that the ground state correlations would decay as power laws, while Yip |^ has found that 
the absorption line-shape would exhibit power-law singularities. 

Motivated by recent experiments Q , we have considered the low-energy excitations 
and ground state properties of finite quasi-ld interacting boson systems. It is important 
to realize that so far the experimental samples are mesoscopic in size |^ and longitudinally 
confined. The experimental consequences of this fact for interacting Id bosons have not been 
fully addressed before. Nevertheless, it is known [^,|ri| that the low-energy spectrum and 
correlation functions of a finite fermionic TLL with open boundaries are very different. In 
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contrast to a system well described with periodic boundary conditions, such as a quantum 
degenerate gas in a tight toroidal trap (a "boson ring"), a bounded fluid cannot sustain 
quantized persistent flows, as a ring would, and its excitations are standing waves. Even more 
important, Eggert et al. |pj]] have shown that the exponents governing correlations in a finite 
bounded fermionic TLL cross over very slowly to the exponents of the infinite system. The 
differences illustrate the effect of boundaries on the quantum critical fluctuations characteristic 
of these systems at very low temperatures, and as we shall argue, they should be taken into 
account when confronting the experiment. Before we proceed any further, it should be noted 
that in this letter the longitudinal confining potential is approximated by two infinite barriers 
at the boundaries. This is not generally the case in current experiments with trapped atomic 
vapors, where the confining potential is usually harmonic. Accounting for this fact within the 
formalism employed in this paper is not a difficult task and does not substantially modify the 
results reported here. Because of space limitations, the necessary modifications to account for 
slowly varying potentials will be presented elsewhere |3C[| . However, it is worth pointing out 
that in recent experiments ^ using microchip traps the confining potential could be tailored 
to a shape very close to that of a square well with very high barriers. Thus we expect that 
our results are directly relevant to future experiments in that type of traps, as well as to the 
"^He systems considered in Refs. [^,|| 

To study the low-energy behavior of a bosonic Id quantum liquid, we rely on Haldane's 
harmonic-fluid approach . How to generalize this approach to deal with boundaries is ex- 
plained below. In general, this method is able to account for the long wave-length properties 
of the system. In this particular case, this means that we cannot describe the properties very 
close to the boundaries, but for most experiments this is not as important as the effect of the 
boundaries on the bulk properties. With these provisos, the low-energy effective Hamiltonian 
as well as asymptotic expressions for the ground state density and one-particle density matrix 
will be derived in this letter. We also analyze finite-size and boundary effects on the momen- 
tum distribution. This quantity is experimentally accessible ||l^,0. Finally, we show that the 
quadratic character of the effective Hamiltonian allows us to derive Bijl-Jastrow wave function 
for the ground state Numerical evidence (e.g. Ref. suggests that our results should 
apply to systems containing even a few tens of particles. 

We begin by assuming that the system is confined by a very anisotropic potential and 
that temperature is much lower than the excitation gap for the transverse degrees of freedom 
(a more detailed analysis of the conditions for (quasi) one-dimensionality can be found in 
|0,||,U,||,|lJ,|l^). A'o spinless bosons in the lowest transverse level of the confining potential 
will be described by the following Hamiltonian, 

H^^J^ dx d,^\x)d^-^[x) + dx£ dx' v{x - x')p{x)p{x'), (1) 

where M is boson mass, p{x) — '^'^{x)'^{x) is the density operator, and ['I'(x), ^^(a;')] = 
5{x — x'), but otherwise commute as corresponds to bosons. The interaction v{x) can be 
totally general as long as it has a short-range repulsive part, and it does not decay slower 
than Confinement in the longitudinal direction, x, is described by imposing that the 

field operator vanishes at the boundaries, i.e. '^{x) — for x — 0, L, where L is the system size. 
As discussed above, apart from this effect, the longitudinal confining potential is neglected in 
what follows. 

For v{x) — g 6{x) the model was exactly diagonalized by Lieb and Liniger Lieb | pO| 
also showed that the low-energy excitation spectrum for g > is adiabatically connected with 
that of a system of impenetrable bosons (the Tonks gas |^, where g — > -|-oo). Girardeau 
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had previously found that the Tonks gas and a system of free spinless fermions have the 
same spectrum. The latter presents a Fermi "surface" consisting of two points, and the 
excitations are particle-hole pairs. The collective modes that exhaust the low-energy spectrum 
correspond to coherent superpositions of particle-hole pairs. Quite generally, a sufhciently 
short-ranged potential will make particles effectively impenetrable when their relative energy 
is small. Therefore, it seems reasonable to expect that the description in terms of the collective 
modes remains good for a fairly large class of models. Yet, different models have different 
high-energy structure and this also leads to sizable corrections to the sound velocity and 
long-distance correlations of the Tonks gas. 

To make the discussion more quantitative, we follow Ref. |^ and work in the density-phase 
representation of the boson field operator, ^'^^(x) = ^ p(x) e~*'^(^\ where [p(a;), e"*"^^^ = 
8(x — x') e~'"^''^\ At low energies the matrix elements of dx4>{x), as well as deviations of 
p{x) from the mean density po (= Nq/L), are small. Thus, if the long wave-length density 
fluctuations are represented as po + n(a;), one has [n(a;), (/((x')] = id(x — x'), i.e. the fields 
n(a;) and (j>{x) are canonically conjugate. The Hamiltonian in Eq. (|l|) can be linearized in 
terms of them to obtain an effective low-energy Hamiltonian, 



-n\x) + -{dM) 

K TT 



^^off = — / dx 



(2) 

where Vg and K should be regarded as phenomenological parameters. They can be determined 
either from an exact solution (when available ^), numerically |13|, or from the 

experiment [^. Once these parameters have been obtained, the above Hamiltonian provides 
a complete description of the low- lying excitations, independently of the details of v{x) 

Next, we introduce another field 6{x) related to H(a;) by dx0{x)/TT = po + ^{x), which 
implies that 9{x) increases by tt every time x surpasses one particle. This means that 

^[6(1) - 6(0)] (3) 

counts the total number of particles in the system. Identifying particle positions with the 
points where 9{x) changes by tt also allows to construct a representation of the full density 
operator |2|, p{x) = dxO{x) Yl^=-oo S{0{x) — mr), which with the help of Poisson's summation 
formula can be written as 

-t-oo 

m— — oo 

The field operator is then given (up to an overall prefactor ||^) by 

+ 00 

^^x) ^[po + n{x)]^ J2 e2™^(") e-^^("). (5) 



m— — oo 



This expression yields an operator that commutes at different points, as can be checked by 
using the mode expansions given below. An anti-commuting operator, which represents a 
Fermi field, can be constructed Q as = 5'^(x)e'^^^-'. The above expressions, Eqs. 

to (|5|), provide us with the tools to compute the spectrum and correlation functions. However, 
we have not yet touched upon the issue of boundary conditions. For open boundary conditions 
(OBC's) we must demand that the field operator vanishes at the boundary, which implies that 

^ e2""«(")=0. (6) 

m— — OO 
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By construction, this is so if 6{0) — 60 n-K-, n being an integer (the precise value of real 
number 9o ^ nn will not be important in what follows); The boundary condition is also 
obeyed for x = L since Eq. (||) relates d{L) to 0(0). The following mode expansions for 0{x) 
and (j)(x) are then obtained (g = rmr/L, for m = 1, 2, 3, . . .) 

9{x) ^ + + (^) ' sin(<zx) [b{q) - hHq)] , (7) 

m = <^o + Yl (-^) ' e-"'/' cos(ga:) [6(g) + h\q)\ , (8) 

which diagonalize the Hamiltonian in Eq. (^): 

iJeff = 5] nu;iq) bHq)biq) + ^{N - N,)\ (9) 

where Lo{q <^ po) = Vg q > 0, and [b{q),b'' {q')] = ^g,g', commuting otherwise. That is, the 
low-energy excitations are linearly dispersing standing "phonons" . The cutoff a ^ p^^ in Eqs. 
(0,H) makes explicit that these expansions are only meaningful as long as we restrict ourselves 
to a low-energy subspace where the phonon wave-length is much longer than p~^. Besides the 
phonons, one can create excitations that change the number of particles. These are described 
by a pair of operators {N,(j)o), which obey [A^, e^'*^"] — e"*"^". In contrast to PBC's, which 
were used in only one pair of these operators is needed because, as mentioned in the 
introduction, one cannot excite quantized persistent flows in a longitudinally conflned system. 
Finally, from Eq. (|^) it follows that the compressibility k = p^'^{dpo/dp) is proportional to 
K. The Tonks gas has K ^ 1 (free spinless fermions), while for weakly interacting bosons 
K +00. 

Using Eqs. (|j7|) we have computed the ground state density for separations from the 
boundaries larger than p^^. To leading order in each harmonic of 27r/9o we obtain 

{p{x)) ^p,n + ^Am , (10) 

I m=l [Pod{2x)\ J 

where d{x) — {2L/tt) |sin(7rx/2L)|, and A^, and 5m are model dependent coefficients HJ]. Since 
in Id the distinction between bosons and fermions is blurred by interactio ns p], the ground 
state density exhibits Friedel oscillations characteristic of Fermi systems [To|,|l2|,|i;|. They are 
induced by the boundaries, which break translational symmetry, and therefore would be absent 
if we had assumed PBC's. For sufficiently strong repulsive interactions [K < 2 P, p5|, |30|] ) , 
these oscillations can be pinned by a periodic potential (e.g. an optical lattice) of wave length 
equal to p~^ . The system would thus undergo a transition to a Mott insulating regime pGf . 

The presence of the boundaries also modifies the structure of correlation functions. We 
have computed the one-particle density matrix at zero temperature. Far from the boundaries 
and for | a; — x' | ;2> ^ , the result reads 



(^'+(X)*(X')) ^ Po 



p-^^d{2x)d{2x') 



d{x + x')d(x — x') 
d{x -f , 



\-oo 

^ ^ ^ —i{m-\-m )TTSgn{x — x )/2 



m.m 



d{x — x') 



2mm K ^TTipo{mx—m'x') 



M2x)r'''M2x')r''''' 



(11) 
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Again Bm,m' , Cm.m' are model dependent (complex) coefficients which cannot be determined 
by this method. Bulk behavior is recovered for |a; — x'| ^ min{a;, }. In 

this limit, the less oscillatory terms in Eq. ( [ll| ) are those where m = m' , which yield the 
expression m Ref. D (27|. Using the same formalism other correlation functions, such as the 
density correlation function, have been obtained. These will reported elsewhere ||30| as only 



the density matrix, Eq.(ll), is needed here 



From the one-particle density matrix the momentum distribution can be obtained : 

n{p^L)^- [ dx j dx' (*t(a;)$(a;')) e*P(^"^'), (12) 
L Jo Jo 

which we have normalized so that (27rpo)~^ / dpn{p) ~ 1. The dominant contribution at small 
momentum is given by the m = m' = term in Eq. ([ll|), henceforth denoted f{x,x',L). In 
the thermodynamic limit L — » oo, this function decays as \x — x'\^^^^^ , which implies that 
n{p <^ Po) ~ with /3 = 1 — {2K)~^ . However, when the system is finite these power 
laws do not hold. This poses a problem to determine the exponent (3 from experimental 
measurements of the momentum distribution. Measuring the exponent is important, not only 
as a means to obtain K and get an idea of the strength of the interactions, but also because [3 
controls the decay of phase correlations: For /3 close to 1 (i.e. X ^ 1) phase correlations will 
decay very slowly, and therefore the system will look like a Bose-Einstein condensate, even 
if, strictly speaking, in Id bosons do not condense. On the other hand, when [3 « 1/2, the 
system will behave as a Tonks gas. 




Fig. 1 - Momentum distribution n{p,L) vs. pL [L is the system's size), for different values of K. It 
has been normalized to its value at p = (Recall that n(0, L) = (poL)^ with /3 = 1 — {2K)^^ . 
Hence n(p, L)/n(0, L) = 7(pL)/7(0), cf. Eq. (p^)). Some results {K — 1, 2, 4) with periodic boundary 
conditions (PBC's) are also shown for comparison. 
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To find the exponent /3, we must resort to a different type of analysis. Let us first note 
that the function f{x,x',L) is homogeneous, i.e. f{x,x',L) — s^^^^ f{sx, sx' , sL). Choosing 
s = 1/L, it follows that f{x,x',L) = L~^/^^ F{x/L, x' /L), for some scaling function 



Introducing this result into Eq. (12) 



n(p<po,L) = {poLf I{pL), 



(13) 



where I{pL) = 2p-l^ d£, d£,' F(C,^') cos[pi(C - £,')]. This function has been plotted in 
Fig. |l| for different values of K. Eq. (|l^) thus implies that the exponent (3 can be extracted 
from a finite-size scaling analysis of the experimental data for n(0, L) collected in systems 
of different size but equal density, po- Girardeau and Wright |l^ have proposed a non- 
destructive method to measure n(p, L). With this technique, one can try to obtain n{p — > 0, L) 
as a function of i, using a blue-detuned laser beam to split a quasi Id trapped atomic gas 
into pieces of smaller size which could be probed independently. Real systems, however, are 
at finite temperature, whereas the previous analysis is for T = 0. At finite temperature, 
correlations decay exponentially with distance Repeating the above analysis n{p, L,T) = 
iPoL)'^r{pL,TL/hv,). Thus, when measuring systems of different size the product TL must 
be kept constant. However, for T < TiVs/imax: where Lmax is the size of the largest system 
considered, quantum fluctuations dominate over thermal fluctuations and the dependence on 
T can be neglected so that our T — analysis should be valid. 

Finally, we derive a Bijl-Jastrow wave function for the ground state of No interacting 
bosons confined in a box with open boundaries. For the Tonks gas, it was pointed out 
by Girardeau that the exact ground state is of Bijl-Jastrow form. Recently, Pham et 
al. ||l^ have remarked that the quadratic form of the effective Hamiltonian, Eq. (||), implies 
that a Bijl-Jastrow wave function is a good approximation to the ground state. Here we 
extend their derivation to OBC's. Expanding H(a;) = Tlo/L + \/2/L J2q>a coa{qx) and 
4>{x) = (f)o + VVLJ2q>o 4'q cos{qx), with [Hg, = iSq^qi and Ho - iVo, Hcs is given by 



2K 



E 

ij>0 



H^ 



qK 



hnVs 
2LK' 



:H^ 



(14) 



where (jiq has been replaced by —iS/dHq, as imphed by its commutator with H^. Thus the 
ground state is a gaussian in momentum space, $o = exp ^—T^J2q>() nq/Sg-fcj, which in real 
space reads 



^o{xi, . . . , xn) = exp < - / dx dx' p{x)IC{x,x')p{x') 



oc 



n 



d{2xi 



n 



d{xi + Xj)d{xi — Xj) 



7r-2L2 



where we have replaced Il{x) by p{x)— po, and then used p{x) = X^i^i 5{x — Xi) and lC{x, x') = 
\og\L~'^T:'^d(x + x')d{x — x')\ / K . The result bears some resemblance to the exact ground 
states of the Tonks gas |Q {K — 1) and the Calogero-Sutherland ||2^ models in a harmonic 
trap. The differences are due to the different confining potential considered. As in those cases, 
however, the wave function vanishes when two particles approach each other, revealing strong 
correlations between the bosons. 



(15) 
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In conclusion, we have extended the harmonic-fluid approach |^ to obtain the low-lying 
excitations and study finite-size and boundary effects in systems of interacting bosons in 
Id that are conflned longitudinally. When analyzing their properties, these effects must be 
taken into account. This has been illustrated by considering the momentum distribution at 
small momentum, which does not behave as the power law predicted for the infinite system. 
Finally, we have also derived a Bijl-Jastrow wave function for the ground state. A more 



complete account will be given elsewhere |30 
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